x=-10:0.1:10;
y=-10:0.1:10;
[X Y]=meshgrid(x, y);
hold on
colormap jet;
subplot(131);
mesh(X, Y, X.^2 + Y.^2);
hold on
contour(X, Y, X.^2 + Y.^2, 30);
subplot(132);
mesh(X, Y, -X.^2 - Y.^2);
hold on
contour(X, Y, -X.^2 - Y.^2, 30);
subplot(133);
mesh(X, Y, -X.^2 + Y.^2);
hold on
contour(X, Y, -X.^2 + Y.^2, 30);
